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Abstract We report results of experiments performed with the KEDR detector at the VEPP-4M e + e~ 
collider. They include precise measurement of the D° and D ± meson masses, determination of the ^(3770) 
resonance parameters, and a search for narrow resonances in e + e~ annihilation at center-of-mass energies 
between 1.85 and 3.1 GeV. 
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1 Introduction 

The paper reports results of three analyses per- 
formed with the KEDR detector at the VEPP-4M 
collider (BINP, Novosibirsk): a measurement of D° 
and D ± meson masses, a measurement of ^(3770) 
mass and widths, and a search for narrow resonances 
in e + e~ annihilation between 1.85 and 3.1 GeV. 

2 VEPP-4M collider and KEDR de- 
tector 

The VEPP-4M collider can operate in the 
wide range of beam energy from 2E= 2 to 12 GeV. 
The peak luminosity in the J/iJ) energy region is 
about 2 x 10 30 cm _2 s _1 . One of the main features of 




C0 2 laser, A=10.59 /im 



Fig. 1. VEPP-4M/KEDR complex with the 
resonant depolarization and the infrared light 
Compton backscattering facilities. 
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the VEPP-4M is its capability to precisely measure 
the beam energy using two techniques Q: resonant 
depolarization and infrared light Compton backscat- 
tering. 

The accuracy of VEPP-4M energy measurement 
with the resonant depolarization reaches 10~ 6 . Such 
measurement requires dedicated calibration runs 
without data taking. The accuracy of the beam en- 
ergy determination for the accumulated data sample 
is dominated by the interpolation between the suc- 
cessive depolarization runs and equals about 6 • 10 -6 
(~10 kcV) in the J/ij) region [2J. 

A new technique developed at the BESSY- 
I and BESSY-II synchrotron radiation sources 0, 
i] was adopted for VEPP-4M in 2005. It em- 
ploys the infrared light Compton backscattering and 
has a worse precision compared to the resonant 
depolarization(50-^70 keV in the J ftp region), but un- 
like the latter it can be used during data taking Q ■ 

The KEDR detector Q includes the vertex detec- 
tor, the drift chamber, the scintillation time-of-flight 
counters, the aerogel Cherenkov counters, the barrel 
liquid krypton calorimeter, the endcap Csl calorime- 
ter, and the muon system embedded in the yoke of a 
superconducting coil generating a field of 0.6 T. The 
detector also includes a high-resolution tagging sys- 
tem for studies of two-photon processes. The on-line 
luminosity measurement is performed with two inde- 
pendent single bremsstrahlung monitors. 




Fig. 2. KEDR detector. 1-vacuum cham- 
ber, 2-vertex detector, 3-drift chamber, 4- 
threshold aerogel counters, 5-ToF-counters, 6- 
liquid krypton calorimeter, 7-superconducting 
coil (0.6 T), 8-magnet yoke, 9-muon tubes, 
10-CsI-calorimeter, H-compensation solenoid, 
12-VEPP-4M quadrupole. 



3 D° and D^ 1 meson masses 

Neutral and charged D mesons are the ground 
states in the family of open charm mesons. Mea- 
surement of their masses provides a mass scale for the 
heavier excited states. In addition, a precise measure- 
ment of the D° meson mass could help to understand 
the nature of the narrow X(3872) state [7- 10 1, which, 
according to some models, is a bound state of D° and 
D*° ll| mesons and has a mass very close to the sum 



of the D° and D*° meson masses. 

Measurement of D meson masses is performed us- 
ing the near-threshold e + e~ — > DD production at the 
-0(3770) resonance with the full reconstruction of one 
of the D mesons. Neutral and charged D mesons 
are reconstructed in the K~tt + and K~tt + tt + final 
states, respectively (charge-conjugate states are im- 
plied throughout this paper). The analysis uses a 
data sample of 0.9 pb -1 . 

The mass of the D meson is calculated as 



\ 



E 2 

beam 



-(x>). 



(1) 



(so-called beam- constrained mass), where -E beam is the 
average energy of colliding beams and Pi are the 
momenta of the D decay products. Since E hcam is 
measured by resonant depolarization technique as de- 
scribed in Sec.[2j its contribution to the D mass error 
is negligible. 

Apart of the M hc , D decays are effectively selected 
by energy difference 



(2) 



where Mi and Pi are the masses and momenta of the 
D decay products. The signal events satisfy a AE ~ 
condition. Therefore, the scale of momenta is tuned 
in order to get (AE) of D decays consistent with 
zero. Furthermore, while calculating M bc , we employ 
a kinematic fit with the AE = constraint. It results 
in a certain improvement of the M bc resolution and 
significantly reduces a contribution of the remaining 
momentum scale error to the measured D mass. 

In order to measure the D mass most efficiently, 
the unbinned maximum likelihood fit procedure is 
used. In the case of D + meson analysis, the likeli- 
hood function depends on two already defined vari- 
ables: M bc and AE. The D° meson analysis em- 
ploys one more variable: the difference of the abso- 
lute values of momenta of D° decay products in the 
CM frame A\p\. We use the fact that the M hc reso- 
lution strongly depends on decay kinematics — it is 
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about three times better if daughter particles move 
transversely to the D° direction (small A|p|), than 
if they move along this direction (large A|p|). Thus, 
the A\p\ variable estimates the M bc resolution on the 
event-by-event basis, improving the overall statistical 
accuracy of the measurement. 

A fit of the event density is performed with D 
mass as one of the parameters in a relatively wide 
region around M hc = M D and AE — (specifically, 
M bc > 1700 MeV, \AE\ < 300 MeV), with the back- 
ground contribution taken into account. The likeli- 
hood function has the form: 

N 

-21og£(a) = -2^1ogp(v i |a) + 2iVlog p{v\a)dv 

i—0 

(3) 

where v = (M hc ,AE,A\p\) are the variables that char- 
acterize one event, p(v\a) is the probability distribu- 
tion function (PDF) of these variables depending on 
the fit parameters a — (M D , (AE),b uda ,b DD ): 

p(v|a) =p sig (v\M D , (AE))+b ud3 p ud3 (v)+b DD p DD (v) . 

(4) 

Here p a i g is the PDF of the signal events which de- 
pends on M D and (AE) (the central value of the AE 
distribution), p uds is the PDF for the background pro- 
cess e + e~ —±qq [q = u,d,s), and p DD is the PDF for 
the background from e + e~ -^DD decays with D de- 
caying to all modes other than the signal one. b uda 
and b DD are relative magnitudes of the background 
terms. The shapes of the p 3 i g , p U ds and p D o distribu- 
tions are obtained from the MC simulation. 

The simulation of signal events is performed with 
the MC generator for e + e" — > DD decays where 
Z?-meson decays are simulated by the JETSET 7.4 
package [jj}. The radiative corrections are taken 
into account in both initial (the RADCOR package 13 1 
with Kuraev-Fadin model .21]), and final states (the 



PHOTOS package |l4(). The ISR corrections use the 
e + e~ — > V(3770) — > DD cross section according to 
a Breit-Wigner amplitude with M = 3771 MeV and 
T = 23 MeV 15]. The backgrounds from the contin- 
uum e + e~ — > qq process and from e + e~ — > DD decays 
are simulated using the JETSET 7 . 4 generator. In the 
latter case the signal D decays are removed from the 
decay table. Full simulation of the KEDR detector is 
performed using the GEANT 3.21 package fltj ]. 

The likelihood fits yield M D o = 1865.05±0.33 MeV 
and M D + = 1869.58 ±0.49 MeV. To obtain the final 
D masses, one has to account for a possible devia- 
tion of the fit parameters M D and (AE) from the 
true D masses and energies. In particular, the cen- 
tral value of M D can be shifted due to the asymmet- 



ric resolution function and radiative corrections. The 
biases are corrected using the MC simulation. The 
final values of the D masses after the corrections are 
M D o = 1865.30±0.33 MeV and M D + = 1869.53±0.49 
MeV, where the errors are statistical only. 





AM D o , MeV 


AM D + , MeV 


Absolute momentum calibration 


0.04 


0.04 


Ionization loss in the material 


0.01 


0.03 


Momentum resolution 


0.13 


0.10 


ISR corrections 


0.16 


0.11 


Signal PDF 


0.07 


0.05 


Continuum background PDF 


0.04 


0.09 


DD background PDF 


0.03 


0.06 


Beam energy calibration 


0.01 


0.01 


Total 


0.23 


0.20 



Table 1. Systematic uncertainties in the D° 
and D + mass measurements. 
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Fig. 3. Comparison of KEDR D meson masses 
with other measurements. The thick and thin 
error bars show the systematic and the total 
errors, respectivel y. T he shaded areas are the 
PDG-2008 values The PDG value for the 
D + is obtained using the measured mass dif- 
ference of the D + and D° mesons. 
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Estimates of systematic uncertainties in the D 
mass measurements are shown in Table [T] Figure [3] 
shows the comparison of our results with those of the 
previous experiments. 

4 -0(3770) mass and widths 



The parameters of the -0(3770) meson were mea- 
sured in the scans of the broad energy region from 
the -0' mass to 3.95 GeV. The analysis is based on 
the luminosity integral of about 2.4 pb _1 . 

Following the similar analyses done previously 18, 
, we parameterize tree-level cross sections of non- 
resonant DD production and 0(3770) resonance as: 



<7™g res (W) = cr° (W ) + a ± (W ) , 
where a ^ (W) — a D £> ■ /3f> ,± , and 



37T 



r.pTfi 



M 2 (PF-M) 2 + r 2 (W0/4' 
where M is the -0(3770) mass, 

T{W) = T tot 



(5) 



(6) 





1 + (H *P D± (W)) 2 


(«o* p r> ( M )) 3 

1 1 / K> „ ... D _ ( 1 / i 1 2 


(Bo*- p D + ( M )) :i 
t- ^— 



r tot is the 0(3770) width, and R is the interaction 
radius. 

The previous experiments [3, have ignored 
the interference between the -0(3770) decay and non- 
resonant DD production. To investigate its impact 
on the extracted parameters of the 0(3770), we pa- 
rameterized the total DD cross section a D £, without 
and with the interference as: 



<Tdd (W) = <t 



i/>(3770) 



{wy 



W (7) 



and 



a DD (W) = \A 



V>(3770) 



Anonres 
' ^-DD 



^i<t>\2 



(8) 



respectively, where \A^ {3770) {W)\ 2 = a i , (3770) (W) 
A"°g res (W) = ^(Ja°g res (W), and is the interference 
phase. 

Parameterization of observable cross section is ob- 
tained by convolution of the a D £> (W) and 0' cross sec- 
tions with the ISR correction function 2l| and beam 
energy spread similarly to Eq. (J9)). 

Fig. H] shows the fits of observable cross section 
without and with the interference. The quality of 
the latter fit is significantly better (addition of one 
fit parameter decreases \ 2 by — 8) and the phase 
is consistent with it. The preliminary results for the 
0(3770) parameters are 
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no interf. x 2 = 17.9/19 d.o.f 
with interf. \ 2 = 9.9/18 d.o.f. 
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> 




3750 



3800 



3850 
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Fig. 4. Fits of observable cross section in the 
0(3770) region: without and with the inter- 
ference. 

• M = 3773.2 ±0.5 ±0.6 MeV, 

• T tot — 23.9 ±2.2 ±1.1 MeV. 

• r ee =294±22±30eV. 

for the fit without interference, and 

• M = 3777.8 ±1.1 ±0.7 MeV, 

• T tot — 28.2±3.1±2.4 MeV. 

• r ee = 312±31±30 cV. 

for the fit with interference. Significant increase of 
0(3770) mass in the latter case is apparent. 
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KEDR(preliminary): 
no interference with interference 
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Fig. 5. Comparison of the 0(3770) mass mea- 
sured by KEDR with other experiments. 

Figure [5] shows comparison of KEDR 0(3770) 
masses with results of other experiments. Mass from 
our non-interference fit is consistent with other mea- 
surements, which ignore the interference, while mass 
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from the interference fit is consistent with the BaBar 
results 20] which also takes the interference into ac- 
count. 

5 Search for narrow resonances 

After the J/ip discovery, search for other narrow 
resonances was performed by several experiments. 
The energy region between J/ip and T mesons was 
explored by MARK-1 at SPEAR &, LENA at 
DORIS 0, and MD-1 at VEPP-4 [24]. The up- 
per limit on the leptonic width of narrow resonances 
obtained in these analyses varies between 15 and 970 
eV depending on energy. The search in the energy 
region below J/ip mass and down to 1.91 GeV was 



performed only in experiments at ADONE 25, 2 



with the upper limit about 500 eV. Recently KEDR 
collaboration has revisited the latter region in view of 
a recent discovery of unexpected exotic states above 
the charm threshold, including the narrow A (3872) 
state, which proved that surprises are still possible 
even at low energies. 

5.1 Experiment description 

The experiment was performed in the beginning 
of 2009 and the results are very preliminary. The 
energy scan started just above J /ip and finished at 
1.85 GeV. The search for narrow resonances was con- 
ducted by automatic decrease of the center-of-mass 
energy by about 2a w (1.4 to 1.9 MeV) steps after 
collection of required integrated luminosity in each 
point. The integrated luminosity per energy point 
varied from 0.3 nb _1 in the upper part of the energy 
range to 0.12 nb _1 in the lower one. The data taken 
at each energy were analyzed in real time. In order 
to improve sensitivity, the integrated luminosity was 
doubled in the energy points with significant excess 
of candidate events. 

The total integrated luminosity of ~ 0.28 pb _1 
was collected. The luminosity for data taking was 
monitored using the process of single bremsstrahlung, 
while the analysis uses the offline measurement based 
on elastic e + e~ scattering in the endcap calorime- 
ter. The beam energy was measured by the Compton 
backscattering technique described in Sec. [2] 

5.2 Data analysis 

The analysis employs three sets of event selection 
criteria with the efficiency to hadronic J/ip decays 
(e h ) varying from 44% to 65%. The e h was measured 
from the observable J/ip cross section and known 
T^l^ ■ Br(J/ip — > hadr) . Below we describe the soft- 



est set of event selection criteria. At the first stage 
we defined track-level criteria to select good charged 
tracks: 

1. Distances of closest approach to the beam in 
the transverse plane and along the beam axis 
are less than 0.5 and 10 cm, respectively; 

2. The deposited energy in the barrel LKr 
calorimeter is above 20 MeV. 

Then an event-level selection, which account a 
calorimeter objects too, was applied: 

1. The total deposited energy in the calorimeter is 
above 400 MeV; 

2. At least two charged tracks in the event satisfy 
first criteria of the track-level selection; 

3. At least one "good" charged track satisfies both 
track-level criteria; 

4. There is a charged track acoplanar to the 
"good" one: |A0 — 7r| >0.15 rad; 

5. Aplanarity of event (sum of momenta trans- 
verse to "event plane") is above 0.1--E 6eQm 

6. There are less than 4 hits in the Muon Cham- 
bers; 

7. lEpf/E^O.5. 

Condition 4 rejects cosmic rays, Bhabhas, and 
dimuon events. Condition 5 rejects radiative Bhabhas 
and dimuons. Condition 6 rejects cosmic ray show- 
ers. Condition 7 rejects two-photon processes and 
hard ISR. 

5.3 Fit results 

The event yield depending on energy is fitted by 
a function that assumes the existence of a resonance 
with mass M R , the leptonic width rf e , and selec- 
tion efficiency e h on top of flat nonresonant back- 
ground. The following parameterization of the ob- 
servable cross section was used: 



a(W) = a + e h J dW dx a R ^ 



hadr(W')- 

W-W 



a w 



(9) 



where 

a R ^ hadr (W) = ^rn.Br(R- 



hadr)-S(W-M R ), 



<7 is the nonresonant background cross section, 
J^(x,W) is the radiative correction function 211 ] . 
G(x) is the Gaussian function. 

The likelihood fits are performed with the M R var- 
ied in 0.1 MeV steps with T^ e -Br(R — > hadr) as a free 
parameter, the energy range in the fit for each M R 
value is M R ±13 MeV. The upper limit obtained with 
this procedure is T* e -Br{R^ hadr) < 100 eV at 90% 
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confidence level. In order to get a more conservative 
upper limit we took into account that: 

• £j/i/>-»hadr/ e e+e-^hadr(3.1 GeV) ~ 1.15, 

• e e + e -^ hQdr (3.1 GeV)/e e+e -^ hadr (1.9 GeV)~1.2, 

• variation of a w could further increase the limit. 

Combining all the factors, we set a limit 

T* e -Br{R->hadr)< 150 eV (90% c.l.) 

in the W range between 1.85 GeV and M J/9 . 

6 Results 

Masses of charged and neutral D mesons are ob- 
tained at KEDR experiment: 

• M D o = 1865.30 ± 0.33 ± 0.23 MeV, 

• M D + = 1869.53 ±0.49 ±0.20 MeV. 

The D° mass value is consistent with the more precise 
measurement of the CLEO collaboration [27], while 



that of the D + mass is presently the most precise 
direct determination. 

The preliminary ^(3770) parameters are obtained 
with standard non-interference parameterization [3, 

• M = 3773.2 ±0.5 ±0.6 MeV, 

• T tot — 23.9 ±2.2 ±1.1 MeV, 

• r ee = 294±22±30 eV. 

The parameters with the interference: 

• M = 3777.8 ±1.1 ±0.7 MeV, 

• T tot = 28.2±3.1±2.4 MeV, 

• r ee = 312±31±30 eV 

Taking the interference into account significantly im- 
proves the fit quality and increases the ^(3770) mass 
by 4.6 MeV. 

A preliminary upper limit for narrow resonances 
in the W range between 1.85 GeV and Mj/ 9 : 

• T*-Br(R^hadr)<m0eV (90 % c.l). 
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